Non-Unitary Spin-Charge Separation in One-Dimensional Fermion Gas 
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In this Letter we report exact results on the infrared asymptotics of one-particle dynamical cor- 
relation function of the gas of impenetrable spin 1/2 fermions at infinitesimal temperature. The 
correlation function shows signs of spin-charge separation with scaling behavior in the charge part 
and exponential decay as a function of the space coordinate in the spin part. Surprisingly, the 
anomalous dimensions in the charge part do not correspond to any unitary conformal field the- 
ory. We find that the fermion spectral weight has a power law divergency at low energy with the 
anomalous exponent —1/2. 
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Significant advance of the theory of strongly correlated 
one dimensional electron systems has beeir made pos- 
sible by two important peculiarities of one dimensional 
physics: universality of low energy properties summa- 
rized by the conformal field theory (CFT) [Jj and the ex- 
istence of exactly solvable strongly coupled models like 
the Hubbard model or the Thirring model 2]. As far as 
their low energy properties are concerned, most physi- 
cally important one dimensional electron systems belong 
to a certain universality class usually referred to as the 
Luttinger liquid (LL) 3j. By bosonization the LL the- 
ory is mapped onto a theory of free massless compactified 
bosonic fields. The compactification radii play the role of 
phenomenological parameters defining the effective the- 
ory. Classification and phenomenology of such effective 
theories is ofteir called g-ology (jj. Due to a relatively 
simple analytic structure of the LL theory, g-ology is a 
powerful calculational tool relating all dynamical infor- 
mation about a strongly correlated system (such as its 
correlation functions, spectral weights, transport coeffi- 
cients, stability with respect to perturbations etc.) to a 
few phenomenological constants, which can be obtained 
from e.g. its thermodynamical properties. For exactly 
solvable models this phenomenological constants are usu- 
ally found from the analysis of thermodynamics by means 
of the thermodynamic Bethe ansatz 



trable spin 1/2 fermions. This model is the infinite U 
limit of the system of spin 1 /2 fermions described by the 
Hamiltonian 
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Despite the success of g-ology, a direct calculatioir of 
infrared asymptotics of correlation functions in exactly 
solvable models has so far been performed only for a few 
special cases 0; IE H IS ^3 and as a general problem 
still remains a challenge. Such direct calculations are 
important for two reasons. First, they provide a test bed 
for the hypotheses underlying the LL theory. Second, 
they give insight into the situations where the low energy 
physics is not described by the LL. 

An interesting example of a physically relevant model, 
where exact calculation of dynamical correlation func- 
tions is possible is a model of non-relativistic impene- 



H = dx 



(1) 

The model can be viewed as a continuum limit of the 
Hubbard model and is the simplest example of interact- 
ing fermions in one dimension. This model has a long 
history. The eigenstates and the spectrum of Hamilto- 
nian ^ for an arbitrary U were found in Thermo- 
dynamics of this model was studied in (l3|. Parameters 
of the low energy LL theory for the Hubbard model were 
calculated in jl,l |. 

The infinite U limit of the model was studied in 
B II II S III III Hi- This is an interesting limit, where 
the ground state of the system becomes infinitely degen- 
erate in spin and the CFT description fails. On the other 
hand, in this limit a specific factorization of Yang's wave 
function occurs 0], opening the way to a direct calcula- 
tion of spectral weights and correlation functions of the 
model through their Fredholm determinant representa- 
tion H . Due to the above mentioned ground state 
degeneracy the limit U oo in correlation functions of 
the model is nontrivial in the sense that it does not com- 
mute with the limit of vanishing temperature T 0. 
Indeed, taking the limit T ^ first one obtains results 
[7], Is], [is' consistent with a naive infinite U limit of the 
corresponding CFT 0,0] . In particular, the one-particle 
momentum distribution function shows scaling behavior 
near the Fermi momentum kp with the anomalous expo- 
nent 1/8 • On the other hand, the one-particle density 
matrix (i.e. the equal time one-particle correlation func- 
tion) decays exponentially as a function of coordinate 
if one takes the U oo limit first. Due to this exponen- 
tial decay the momentum distribution function shows no 
scaling near kp. 
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In this Letter we report our results on the asymptotics 
of the one-particle dynamical correlation function 



G{x,t) = U\{x,t)il;^iO,0) 



(2) 



when the limit J7 ^ cxd is taken before the limit T — > 0. 
Due to the rotational invariance of Hamiltonian the 
correlation function of spin down fermions coincides with 

®- . . . . 

We begin our discussion with first recalling how the 

infinite U limit in correlation functions is taken at zero 
temperature and then we discuss the opposite order of 
limits qualitatively. Then we present our exact result 
for the asymptotics of the correlation function Q and 
discuss its structure and physical meaning. 

Assume that in the ground state of the system there 
is a finite fermion density Uc — n-\ + n^, controlled by 
the chemical potential ^. At any finite U the low energy 
physics of this model can be described by the LL the- 
ory. According to this theory the low energy spectrum is 
described by the effective Hamiltonian 



Hr + Hs 



(3) 



where He describes charge fluctuations, describes 
spin fluctuations and [Hc,Hs] — 0. Independent dynam- 
ics of spin and charge is called spin-charge separation 
0. Hamiltonians He and Hg are, in fact, Tomonaga- 
Luttinger Hamiltonians corresponding to c = 1 confor- 
mal field theories with global U{1) symmetry. Hamilto- 
nian Hg also possess global SU{2) symmetry inherited 
from the rotational symmetry of This additional 

constraint on Hs fixes the structure of its spectrum com- 
pletely up to a non- universal constant v^, which is the 
spin propagation velocity depending on the microscopic 
parameters. The fermion operator in the LL theory is 
represented as a sum of anticommuting "left" and "right" 
fermion fields 



Akpx 



=T,i (4) 



where kp is the Fermi momentum. The "left" and "right" 
fermion fields are products of spin and charge vertex op- 
erators 



(5) 



where the spin operators satisfy S'L(fl).| = '^l{R) y ^P" 
erators S and C commute with each other and satisfy 
[S, He] = [C, Hg] = 0. Formulas ©, (B and © only ap- 
ply to the low energy properties of the system. That is, 
the energy must be much smaller than a certain cutoff 
scale A, determined by /i and U. 

At zero temperature, the correlation function ^ pre- 
dicted by the Luttinger theory has the form 

G{X, t) = §(2; - Vst)e{x - Vet, X + Vet)e-'^^'' 

+2>{x + Vst)Q{x+Vet,X-Vet)e^^'''' (6) 



where 

Kx-Vst) = l^S^^{x,t)SR{0,Q)) = 
and 



const 



{x-vsty/^ 



(7) 



C(a; — Vet, X + Vet) 



Cl{x,t)CR{Q,Q) 
const 



{X - Vet)'^^-{X + Vet) 



2Aa ' 



(8) 



where Vg and Ve are the propagation velocities of spin 
and charge excitations respectively. Correlation func- 
tions (|7|), (IHl do not contain any length parameter. This 
refiects the scaling invariance of the LL theory. Accord- 
ing to O the spin operator Sr^^ is a right mover with 
fixed anomalous dimension 1/4. This is a consequence of 
the constraint imposed on the effective Hamiltonian Hs 
by the global SU{2) symmetry. At the same time, the 
correlation function of the "right" charge operator Cr^ 
splits into a product of left-moving (depending on x+Vet) 
and right-moving (depending on a; — Vct) parts with dif- 
ferent positive anomalous dimensions A and A. These 
dimensions are determined by the Luttinger parameter 
K, which can be found from the charge compressibility 
K = due/dfi of the system as follows 



One has Q 



Ae 

A, 



16 V Vk 
16 V Vk 



(9) 



(10) 



For any finite U parameters Ve and k can be calculated 
from the thermodynamic Bethe ansatz. For infinite U 
fermions become impenetrable, which from the point of 
view of thermodynamics means that they satisfy the 
Pauli principle independently of their spin orientation. 
For such a system parameters Ve and k are calculated in 
the same way as for a system of spinless non-interacting 
fermions |l6L llTj . Calculating the right hand side of Q 
for spinless fermions one finds 



K 



(11) 



For the anomalous dimensions in the charge sector this 
gives 
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32' 



Ae 



1 

32' 



(12) 



Correlation function for impenetrable fermions in the 
form with anomalous exponents in the charge sector 
given by (|12|l was suggested in Q . It was also discussed 



3 



in [isj , where the U oo hmit was analyzed microscop- 
ically. 

Next, consider the infinite U limit taken at finite tem- 
perature. For large U there happens a separation of en- 
ergy scales in the problem. This separation is controlled 
by a small dimensionless parameter kp/U. While /i re- 
mains the characteristic cutoff energy in the charge sec- 
tor, the characteristic cutoff energy in the spin sector 
becomes of the order of iikp/U. The same happens with 
the velocities. While the propagation velocity for the 
charge Vc — > as U — > oo, the velocity of spin exci- 
tations vanishes, Vs — > 0. In the limit of infinite U the 
spin degrees of freedom completely loose dynamics and 
the ground state of the system becomes infinitely degen- 
erate with respect to spin flips. From the point of view 
of the low energy physics it is interesting to consider a 
situation where 



(13) 



In this situation all spin configurations have equal sta- 
tistical weight and the dynamics of the system is effec- 
tively averaged over spin configurations. At the same 
time, the charge degrees of freedom are not strongly af- 
fected by temperature, because the latter is much lower 
than the corresponding cutoff energy /i. Thus in the tem- 
perature range \i'6\ correlation function (|2J) should not 
strongly depend on temperature and its behavior should 
be in this sense universal. The infinite U limit of the 
correlation function taken at a fixed small temperature 
T ^ fi thus corresponds to taking the infinite U limit 
first and then taking the limit T —^ 0. One can try to 
conjecture the structure of the one-particle correlation 
function ^ in this situation using formulas (0, JHJ and 
the general results of the CFT. For finite temperatures 
and for x > Vs/T, the spin part of the correlation func- 
tion predicted by the LL theory is 



%r{x - Vst) 



(14) 



This formula is valid for temperatures smaller than A^. — 
fj,kp/U ~ kpVa, that is when linear description of spin 
waves is appropriate. For |a;±Wci| < Vc/T the charge part 
of correlation function is given by ||SJ). For T > Kg all spin 
degrees of freedom become saturated and the spin part 
of the correlation function should become temperature 
independent. Thus, replacing the temperature T in the 
right hand side of (|14|l by the crossover scale A^ we get a 
prediction for the correlation function at T > A^. Making 
this substitution and taking the limit of vanishing Vs we 
get for the fermion correlation function 

G{x, t) = ${x)t{x + Vet, X - Vct)e-'^'''-' 

+%{x)t{x-Vct,x + Vct)e'^'"', (15) 

where function C is given by (|SJ) and 



Here 7 is a dimensionless constant. The problem of this 
picture is that it assumes spin-charge separation in the 
LL sense even though the spin degrees are strongly ex- 
cited. To see that this is not the case, compare the cor- 
relation function of spin up fermions calculated in two 
different ground states. Let one ground state be the infi- 
nite U limit of the ground state of Hamiltonian and 
another be the fully polarized spin up ground state. Both 
states have the same statistical weight and equally con- 
tribute to the correlation function. In the first state the 
anomalous dimensions of charge part of fermion corre- 
lation function are given by (|12() . In the second state, 
due to the Pauli principle fermions do not feel the inter- 
action and the anomalous dimensions are those of non- 
interacting fermions, that is A = 1/2 and A = 0. It is 
not, therefore, clear whether after averaging over all spin 
configurations the anomalous dimensions in (|15|l should 
coincide with (|12|l or, even, whether the correlation func- 
tion should have the structure (|15(l at all. 

In paper il4i | it was shown that the correlation function 
(O can be expressed in terms of the Fredholm determi- 
nant of a linear integral operator V via 



G{x,t) 



The kernel 



87ri 



dz 

—F{z)B__{z)Aet{I + V){z). 



(17) 



k — p 

of the operator V is defined on the square [—kp,kp] x 
[—kp,kp], where kp = ^fji. Functions entering p8|l are 
defined as follows (z = e"*) 



a/tt 



(fe)/2 



(19) 



00 g~r(p) 



Tr(j) — k) ' 



e+(fc) = \ [(l-cos77)e"W^o(fc)+sin7?], (20) 

Eo{k) = p.v. / dp 
T{k) ~ ik^t — ikx. 
Function B__{z) is 

B—iz) = /" ^ dke-{k){i + V)-^e-{k) (23) 



(21) 
(22) 



and 



F{z) = 1 



2- z 2z - 1 



(24) 



= const X e 



(16) 



We performed asymptotic analysis of Eq. fT7| 
by asymptotically solving the corresponding Riemann- 
Hilbert problem by the techniques described in [l8| . 
Detailed calculations are involved and will be presented 
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elsewhere Our main results for the correlation func- 
tion ^ are as follows. For x,t ^ +00 and x/t = const 7^ 
we have 

G{x,t) 



{x-2kFtf^{x + 2kFtf^ 

'^^ — kp-x In 2/7T ^—i{kFX — 4>— ) 

~ {x + 2kFty'^{x-2kFt)^^' 



(25) 



Here is a constant, which is calculated explicitly in [if 
The phases 4>± are given by 



= -2Im 



InT 



iln2 
2tt 



In 2 

ln(2fci.a;±4fc^t) (26) 

TT 



and the anomalous dimensions are given by 



2 8 



In 2 

TT 



This result (|25|l is consistent with the previously calcu- 
lated equal time correlation function . The case x = 
is special. In this case 



G(0,<) = 



77' 



(28) 



^t\n{klt) 

where S' is a constant. 

The form of the correlation function H25|l is essentially 
the same as of (|15|) and shows clear signs of spin-charge 
separation. The exponentially decaying part can be at- 
tributed to the non-propagating spin mode. The alge- 
braically decaying part of the correlation function only 
depends on combinations x — 2kFt and x + 2kFt which 
correspond to propagation with velocity 2k f- This is ex- 
actly the velocity of the charge propagation in the model. 

Remarkably, the anomalous dimensions given by H27|l 
are completely different from dimensions H12|l . obtained 
in the infinite U limit of the LL theory. While the form 
of the correlation function (|25|l suggests that the fermion 
operator can be represented as a product of a spin oper- 
ator and a charge operator, the emergence of a negative 
scaling dimension A in H27|l indicates that the conformal 
field theory describing the charge sector should be non- 
unitary. Answering the question of whether such a CFT 
exists is a subject of further investigation. In particular, 
it would be interesting to calculate the three point corre- 
lation functions involving two fermion operators and the 
charge density or the spin density operator. 

Note that the phases 0± (|2t)|) instead of being constant 
are logarithmic functions of the light cone coordinates 
X ± Vet. A similar phenomenon occurs in the theory of 
classical non-linear wave equations [ioj. 

While the exponential decay of the correlation func- 
tion in X direction smears the power law scaling of the 
momentum distribution function near the Fermi points, 
the power law scaling can be observed in the tunnelling 
density of states, which is given by 



where tilde denotes the Fourier transform. From our re- 
sult (|2Hl we find 



A{co) 



(30) 



The divergence of the tunnelling den sity of states at zero 
energy is similar to the one found in [l^ but the numeri- 
cal value of the tunnelling exponent 1/2 is different from 
3/8 calculated in [l5| . 
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